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ABSTRACT 

An accurate method of parameter estimation for the mathemati- 
cal modeling of semiconductors using the Ebers-Moll equations is 
presented. Its usefulness is apparent in estimating parameters 
to be used in computer circuit-analysis programs that have been 
developed. The Ebers-Moll models were modified to better 
represent the actual characteristics. The least-square-error 
methods presented for estimating the parameters are easy to pro- 
gram to the digital computer and result in parameters that are 
quite accurate in describing the actual characteristics. This 
procedure yields solutions for parameters that are quite helpful 
to the engineer in solving electrical circuits involving p-n 


diodes and junction transistors. 
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IT. INTRODUCTION 


In this age of computers, several excellent computer circuit- 
analysis programs have been developed which solve the equations 
describing electrical networks with remarkable accuracy. In order 
for the computed results to agree with actual circuit performance, 
the procedure of device modeling must furnish accurate models 
which are computationally efficient. As stated in a paper 
published by Design Automation, Inc. ale "For all its advantages, 
the computer analysis of a circuit can be no better than the data 
that are fed into it. Poor models or improperly assigned values 
yield inaccurate results. In fact, once a circuit-analysis pro- 
gram iS running on a computer and the designer has developed the 
skill to use it, the choice of semiconductor models and procedures 
for finding numerical values for the models' parameters are the 
main stumbling blocks to applying the program. The parameters of 
a real switching transistor or real diode vary as functions of 
voltage, current and temperature. Modeling is most accurate when 
the device characteristics are measured under conditions of 
voltage, current and junction temperature similar to those under 
which the device is expected to operate." After a model has been 
chosen, it is essential that the various parameters must be 
computed easily and accurately. With nonlinear devices, such 
as diodes and transistors, the model eyuations must give a 
mathematical description of the curve so that for a given voltage 


the current is known or vice versa. For this analysis, the 


Ebers-Moll [ 2] models (modified slightly) for the diode and the 
transistor are usually employed. Using these models, methods 

are needed to determine the various parameters necessary for 
computation using any of the several computer circuit-analysis 
programs (TRAC, CIRCUS, NET I, NET II, and SCAN). These models 
Cover both small-signal linear and large-signal nonlinear 
operation. The analysis in this thesis is limited to determining 
model parameters using steady-state characteristics. 

In May, 1965, engineers at North American Aviation, Inc. 
completed analysis on the least-square-error curve fit for the 
semiconductor diode model parameters for the SCAN and TRAC 
computer circuit analysis programs. The modified Ebers-Moll 
model for the diode was presented and the necessary assumptions 
were included. 

Design Automation, Inc., in June 1967, proposed a modeling 
procedure for transistors and diodes for computer-aided nonlinear 
Circus t analysis ‘1 . The models presented were essentially 
the Ebers-Moll models with additions. The analysis by Design 
Automations, Inc. was a for the NET I Network Analysis 
Program. Several improvements were made to the existing models. 

In June 1968, the junction transistor model (modified Ebers- 
Moll model) was presented by E. Steele for the TRAC computer 
analysis program er In March 1969 engineers at North American 
Rockwell Corporation presented, in an internal letter, a proposed 
method, using least-square estimation, of deriving the various 
parameters by analyzing the two junctions of the transistor 


separately. 
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In this thesis the problem of estimating the parameters of 
the Ebers-Moll models is considered. Questions investigated 
are: (1) determination of which regions of a semiconductor device 
are most useful for estimating each parameter of the model, 

(2) development of computer techniques of parameter estimation 
using least-square-error and normalized least square error 
solutions, and (3) estimation of the accuracy of the estimations. 
The ultimate question is whether the proposed models give an 
accurate representation of the actual performance of the semi- 
conductor devices. 

In Chapter TI the. modified Ebers-Moll models for the diode 
and the transistor are presented. These models allow accurate 
calculations to be made of the semiconductor device parameters 
necessary for the computer circuit-analysis programs. 

In Chapter III the two methods of least-square-error 
solution are presented. One method is based upon the usual 
least-square-error solution and the other is based upon a normal- 
ized least-square-error solution as developed by Werther and 
parker [14 ]. 

Procedures for calculating the semiconductor device parameters 
are presented for both the diode and the junction transistor. A 
discussion of the useful regions of parameter determination and 
the typical values expected from this analysis is included. A 
method of accumulating measured data and correcting estimated 
values accordingly in also discussed. 

The parameters for typical diodes are calculated in the 


Appendix, using the two least-square-error methods. A comparison 


ale: 


of the results of these two methods is made and includes 
calculation of the error aiter the model parameters have been 
estimated. 

The computer programs used for the solution of the model 
parameters are given following the Appendix. Recursive equations 
for both least-square-error methods are used. The FORTRAN IV 
programming language is used for the computer programs. Most of 
the subroutines were developed for this analysis, but several 
were existing subroutines in the Naval Postgraduate School 


Computer Facility subroutine library. 
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II. EBERS-MOLL MODELS 


Ebers and Moll described the diode model and the junction 
transistor model, under DC conditions, using the small-signal 
parameters to relate to large-signal behavior [2]. These models 
described all regions of operation. Most of the other proposed 
models are mathematically equivalent and generally less accurate 
than the Ebers-Moll model [s]. Both the diode and the transistor 
models are considered below. The Ebers-Moll equations are modi- 
fied by considering a term called the emission constant M for the 
G@tode dnd Mp and Mc for the transistor. Also linear resistanees 
are added to account for the deviations from the ideal models. 
The Ebers-Moll models relate the semiconductor device character- 
istics in terms of several parameters, and generally provide very 
Jittle insight to the actual physical processes within the devices 


represented by these models. 


A. DIODE P-N JUNCTION MODEL 
For an ideal p-n junction diode, the current and voltage 
have the following relationship 
T = 1, : [exp @/kt) "aaalll (2. 

where I. is the reverse saturation current 

V is the voltage across the terminals 

q is the charge on an electron 

k is the Boltzman constant 

T is the absolute temperature (room temperature-298°K) 


kT/q equals 26mV at room temperature. 
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In practice very few p-n junctions can be described using the 
ideal diode relationship. In order to account for deviations 
between the theoretical and actual characteristics the ideal p-n- 
junction diode needs to be modified by placing linear elements in 
parallel with or in series to the ideal diode. A large parallel 
resistance called the reverse diode leakage resistance is used to 
represent the slope of the reverse characteristics, and a ste 
series resistance called the diode ohmic series resistance is used 
to represent part of the voltage drop across the p-n junction. 
Also, a modification to the Ebers-Moll model is the quantity 
M called the diode emission constant, with a value between 1] and 3, 
which changes the exponential dependence to more exactly describe 
the actual behavior, because of the non-ideal operation [3]. 

The wodiiied ideal diode equation becomes 
I = I, jexp(qV/MkT) - ii (2.2) 

The DC diode circuit is given in Figure ] using the above model. 
The various quantities included in Figure ] are defined as 

Vj, Ij - the measured diode terminal voltage and current 

Ree - tHe combined body, lead and contact resistance 

Ry - the reverse surface-leakage resistance across the junction 

D - the modified ideal diode 

I - the ideal diode current given in equation 2.2. 
From the measured data of current and voltage, the parameters 
I_, R_, R, and M are to be determined to describe a specific 


S S 1 


diode or p-n junction. 
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FIGURE 1 


DC P-N Diode Circuit Model 


LS 


B. JUNCTION TRANSISTOR MODEL 
For the ideal junction transistor, Ebers and Moll proposed 
the following current-voltage relationship assuming a pnp 


transistor 2 |. 





Ip = +EO xp (V,./kT) -1 - xp (GV./kT) -1 (2.3) 
N 


| Il 
{exp (qVp/kT) -1} + era Xp (qVc/kT) -l (221) 
N 





(2%) 
where I, = the emitter current 
Ta ethewco) Lact om semicre nt 
Ve =@ehe emitter-to-base voltage 
Vc = the collector-to-base voltage 
Nien cle common-base normal-mode DC current gain 
J = the common-base inverted-mode DC current gain 
Ipg = the Sdluration ¢yrrent of the emitter junction wii 
zero collector current. 
leg = the, saturation current of the collector junction wile 
zero emitter current. 


From this the following is true [3] 


Teg = Teop/C- Ay YX) (2.6) 
Teg = T¢9/-AX yy) (259) 

where Inc = the emitter-to-base p-n junction saturation current 
Tec = the collector-to-base p-n junction saturation current. 
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In order to make the exponential dependence more exact, two 
Gudntities Mp and Mc are needed similar to the diode exponential 
emission constant. The quantity Mp is the emission constant for 
the emitter-base p-n junction and the quantity Mc is the emission 
constant for the collector-base p-n junction. The ideal transis- 


tor equations then become 


Ip = Tgg [exp (@Vp/MekT) -1|- XrI¢g [exp (aV¢/MckT) -1] (2.8) 
and 
Io = Igg [exp (aVe/MckD) -1}- Xgl eg [exp @Ve/MgkT) -1] - (2.9) 


Linear circuit elements are then added to modify the ideal transis- 
tor to an equivalent element to describe the actual relationship, 
Similar to the diode formulation. The transistor equivalent 
circuit is presented in Figure 2 [2 | . The elements of Figure 2 
are defined below: 


fe ip + la@the base cumrenc 

I5 = the current generator across the collector-base junction 
= las exp (qV./M, kT) - - ee [exp (qV,./M,kT) -1 | 

ime — ene Current generator across thesenitter-base junerion 
= Ipg [exp (GV—p/MekT) -]]- %Icg [exp (GVe/MckT) -1] 

Rep = the base spreading, bulk and contact resistance 

R = the emitter bulk and contact resistance 

Rec = the collector bulk and contact resistance 

Re — ehe emitter-base junction ohmic Leakasevrecsceaner 


Ro = the collector-base junction ohmic leakage resistance 


ie 


=> (ent 
a 
Re Re R mE 
~— — 
| x 
k “ Ne 
NN, Ree 
ie 
8 
PEGUKE 2 


DC Junction Transistor Circuit Model 
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The equivalent circuit for a pnp junction transistor may be used 
for a npn transistor by changing the polarities of Vc and Vp and 
wee directions of 1; and I>. 

The specification sheets for transistors usually specify 


common-emitter DC current gains H and H rather than the 
FE, FEs 


common-base DC current gains Xn and %,. The relationship 


between these &'s and Hpp's are 
ro N= Hee / Wars (27008 


and 


x 


I! 


1 = Hpe./(Hpp +) - (2.11) 
Al: 1 


A simplified ideal model of the junction transistor is given 
iieereure 3, Were’ 1, and I> are replaced witthgdiodes and 
dependent current sources. The resistances Rp and Roe are assumed 
to be very large and Rpp, Rec, and Ref are assumed to be small. 
For Figure 3 the elements are defined below 


Ipp = Ipglexp (qVg/Mgk7) -1] (2.12) 


and 


Icr = Icg [exp (qVc/Mckt) -1] (2.13) 
From measured data, the parametersRcc, Ref, Rpp> Me. Mc, 
fee. Icc, BEES and Hee. are to be determined to describe a 
specific transistor. 
The complete model including the AC components for the diode 
is given in Figure 4. The AC components added to the diode model 


are 


C, - the junction transition capacitance 
and 


Se 


Lee Ler 





FIGURE 3 


Simplified Transistor Model 
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ie 
= 


FIGURE 4 


Complete Diode Model 


el 


Cq - the junction diffusion capacitance. 

In Figure 5 the AC components have been added to the 
transistor model so all operating conditions except the break- 
down regions for the junctions can be accounted in the model. 
The added components are 

Cip - the emitter-base transition capacitance, 
Ci+c - the collector-base transition capacitance, 


Cap - the emitter-base diffusion capacitance, 


Cac =sime colleetor=base ditfus#en edpaci tance’. 


The transition capacitances are functions of the junction voltages 
and the diffusion capacitances are functions of the emitted 


currents from the emitter and from the collector [1]. 


ad 


> <— 
ed es 
ue 
‘¥ P 
WV 
Ree Re Re 
: Cac Cd E 
AL = a 
Cc 
\ 
R BB 
\ I 
FIGURE 5 


Complete Transistor Model 
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III. LEAST-SQUARE-ERROR SOLUTIONS 


The method chosen to compute the various parameters of the 
models presented is the least-square-error estimation. This 
solution is for a set of overdetermined simultaneous linear 
equations. There are two least-square-error methods presented. 
The first method is the usual least-square-error solution 
utilizing the concepts of the pseudo inverse as developed by 
Penrose [6] ana presented by Werther [ u]. The second method is 
the normalized least-square-error solution as developed by 
Werther and Parker [u]. A recursive equation is presented for 
each method to facilitate the use of additional data points in 
the calculation. These two methods and the estimations are 
compared with measurements to determine the usefulness of each 
method. 

The material for the least-square-error estimation and the 


normalized least-square-error estimation was obtained from 
Werther [4 | , 


A. LEAST-SQUARE-ERROR SOLUTION 

The least-square-error solution is a very common method of 
solving a set of m simultaneous linear equations in n unknown, 
m being greater than n. This is the overdetermined case for 
linear equations which is common when measurements are made to 


find a small number of unknowns. The least-square-error solution 
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solves the Pelantoneneee 
Ax = b (3.1) 
where A is the m by n matrix of coefficients, 
Deis she mebyel avector of constants, 
and 
x is the n by 1 vector of unknowns. 
A unique solution for x is obtained by using the concept 
of the Penrose pseudo inverse. This yields the exact solution 
if it exists, or the least-square-error approximation if the 
exact solution does not exist. The Penrose pseudo inverse is a 
n by n matrix and is identified as At. The least-square-error 
solution is 
= aty. (3.2) 


If the matrix A has rank n the least-square -error solution is 


zal 
< =[aT “aT (3.3) 


-1 
where A! is the matrix transpose of the matrix A and [ata] de 


the matrix inverse of the product of Ne and A. [In order for the 
rank of A to be n, the rows of A must be independent. This 
solution is easily obtained on the digital computer with available 
matrix operation subroutines. 
For the recursive relationship, the least-square-error 
solution for the first k points is 
Me = Bey 2y (3.4) 


where Z, is equal to by, which is the vector of k points, 


lThe bar under a lower-case letter represents a column matrix 
or vector. 
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Ay is the k by nm coefficient matrix 
ale 
P, is theminverse ot the mera xe, 
and 


Q is the least-square-error solution. 


Now to consider the k+l equation, ae 


to the A matrix and Z,,) is the new data for the vector b,j}. 


Therefore 
A 
= k 
rT = 
Ab 
a 
and 
Zz 
he ball 
bee = 
2k+1I] . 


The k+l least-square-error solution is 


A 
b T 
Xe = Pace Zt 


where P,.4)] is the inverse of the matrix | Ay+yAq41 | - 


The resultant recursive equations are 


N /\ Pra 
Xt] = Xk + A kt -a xn) 
l+a Pra 
and 
p -p _ Pka a Pk 
kt+1 ae 
l+a*Pka 


B. NORMALIZED LEAST-SQUARE-ERROR SOLUTION 


The normalized least-square-error solution is a weighted 


is a row of new coefficients 


(3.5) 


(3 . 6) 


(3.7) 


(3.8) 


(3.9) 


least-square-error solution with weighting factors chosen such 
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that the solution lies as close as possible to all geometric loci 
described by equation 3.1. 
The matrix of weighting or normalizing factors is a diagonal 


m by m matrix called W. The factors w;; are given by 


Nn 
h 
Wee > 2332) Se logs: Sin (3.10) 
il 


This weighting matrix is multiplied times the matrix A and 


vector b resulting in two new quantities 

A* = WA (3.521) 
and 

i ee or, 12) 


The normalized least-square-error solution is then 


th 
Qe = [axTax] “aeTpe. (3.13) 


The recursive relationships for the normalized least-square- 


error solution are 


A ay re are a* 


NneAS mee" (3.14) 
ANk oer l+a*Pf ax CAT 7 axe.) 
and 
Pi axa p# 
oe =p - — a : Go) 
Using the fact that 
L 
3 a af —% 3,6 
a) eee _— 


ag 


and 


at = (a’a) a (3.17) 


the recursive relationships become 


fe 
rx A P, a 3.18 
X*k+] = X*,4 So ee ~ alxx ) (3.18) 
pare a ioe 
and 
ke T 
inten k (3.19) 
* = Ps = ° . 
k+1 7% 
la le a 


The recursive relationships for both the least-square-error 
solution and the normalized least-square-error solution are 
easily adapted to the digital computer and a sample computer pro- 


gram for each method is found after the Appendix. 
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IV. DETERMINATION OF SEMICONDUCTOR PARAMETERS 


Using the models presented in Chapter II, the model parameters 
are calculated utilizing the least-square-error method. Param- 
eter calculation is accomplished using the above equations and 
the equivalent circuits. In order to use the least-square-error 
methods, the nonlinear equations characteristic of these semi- 
conductor devices must be converted to a set of linear equations, 
which involves several assumptions. The solution is made on the 
nonlinear equations where these equations can be assumed to be 
linear. The error associated with these assumptions is so small 
it is considered negligible. 

These parameters enable the user to employ these models 
in existing computer circuit-analysis programs and obtain 
computed results that very accurately represent the actual per- 
formance. The data for parameter estimation may be obtained 
using the measured current and voltage values or values taken 
from device specification sheets. When specification sheets 
are used, the values of the voltage should be read to several 
millivolts in order to obtain accurate results. 

Since the data involves experimental error and is statistical 
in nature, as many data points as possible should be used in the 
estimation of the semiconductor parameters. Since the measure- 
ments are independent, in general, the rank of the matrix A in 
equation 3.1 is equal to n, the number of parameters. Since the 


overdetermined case must have the matrix A with a rank of at 


least n, the solution presented in Chapter III may be used. 


ey 


A. LINEAR EQUATION FOR DIODE PARAMETER ESTIMATION 
The model for the semiconductor diode is presented in 
Figure 1. In order to determine the values of the parameters, 
the model equation has to be simplified and made linear to be 
compatible with the least-square-error methods. Using Figure l, 
the following can be stated: 
Le = ee (4.1) 


Or 


I, = I, [exp (QV/MkT) -1j + V/R,- (4.2) 


The tirsterestricticn wis, thai) (the .measured ste midis: 
voltage) be positive (Vj 70) . With this restriction, 
exp (qV/MkT) >> 1 because even with small positive measured 
voltages (for example 0.1 volt) and with q/MkT = 0.025, the term 
exp (qV/MkT) =54.99. This is much greater than 1. Since R] is so 
very large compared to V, I >> V/R,, and equation 4.2 is simpli- 
Pied 0. 


I, = I,exp (qV/MkT) (4. 3) 


which by rearranging becomes 
V = MkT/q In(15/T,) . (4.4) 
Referring back to the model 


- = + “ 2 
Naw + Lei (4.5) 


and substituting 4.4 into 4.5 results in 


\- 


4 (4.6) 


ACME Pa ae 
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This equation is simplified by letting Ag = MkT/q and 


a - Wi_., obtaining 
V5 = Ag ele) + Ao in(I5) 33 I4R, (4.8) 
Vi = A3 + Ap 1n(T5) + I4R,, (4.9) 


perc Ase= Aj In(A)). This approach is currently used by 
North American Rockwell Corp., among others. 

In the above equation, V5 and I; are the measured forward- 
voltage and current data points for the diode. The other three 


quantities are constants to be determined using the least-square- 


error procedure. For the least-square-error method 





oT 
A = (4.10) 
Tn 
b= (4.11) 
and 
A3 
2 = a (4.12) 


where & iomime least-sGldmeserren colmiaion. Aiter A>. SAgmand Kc 
are determined, the parameters of primary interest are calculated 


using the following relationships: 


ou 


M = Agg/kT, @.13} 
A, = exp(A3/AQ) . (4.14) 

and 
I, = 1/A,. (4.15) 


The three parameters M, Ig and Rg are then calculated and 
can be used in the circuit-analysis computer programs. 

For determination of R], (the reverse diode leakage resist- 
ance due to the imperfections, surface effects, and thermally 
generated carriers in the depletion region) the same equation is 
used as the one for the determination of the other diode 
parameters. For this case, data points are used from the reverse 
characteristics of the diode; thus Vj and Ij are both negative. 
The basic equation is 4.2. By making an assumption that 
exp (qV/MkT) <1 for V less than -0.2 volts,* equation 4.2 is 


reduced to 


— 
II 


7 “la + V/R} - (4.16) 


solving for V, 


V = (13 + Is)R (4.17) 


and substituting this into equation 4.5 results in 


\. 


; dhe + Tae, + I5R,. (4.18) 


Rearranging, 


<= 
| 


eer (pe aseere TORN (4.19) 


ewith V = -0.2 volts and MkT/q = 0.026, exp(qV/MkT) = 4.5x107> 


and most measured data is less than -0.2 volts. 


ie 


couemiemany | Ro) > aewancere Rk), = AE results in 


V; = IyAs + Ag. (4.20) 


This is solved using the least-square-error procedures with 





T 
A (4.21) 
i 
b (4.22) 
x = ee (4.23) 
x* |e 


In the above procedure Vj; and I; are the measured voltage 
and current in the reverse direction of the diode. From the 
results, R] isadetermimed after an assumption isemade. Exam) the 
assumption that R] >>Rs, R] is then equal to As since As = R1 + Rg 
and i 3 is so small it can be neglected. A value of reverse 
seriration current, I., may be foumd from th#s data ‘but this 
value is not a good value for this current, especially for silicon 
diodes. For a 1IN540 diode, the reverse saturation current is 
determimed from the reverse characteristicswalong with@Ry] -- The 
value obtained is then placedas a constant into the equation used 
to derive the forward parameters. The other two parameters are 
derived to see if this reverse saturation current resulted in 
parameters in the useable range. It is found that the value of 


Rs is -1.46 ohms and the value of M is 3.67. These values are not 


53 


realistic numbers to be used in computer circuit-analysis pro- 
grams. These values are determined in the Appendix. The reverse 
characteristics of silicon diodes tend to be more difficult to 
model, since the ideal diode law breaks down completely as a 
description for the reverse silicon-diode Ponay ior. male In most 
p-n junctions the reverse current is not completely voltage- 
independent, due to the surface leakage. R, is the only reliable 
parameter derived from the reverse data. 

Several regions of the forward characteristics were investi- 
gated to discover where the best parameter estimations can be 
found. The standard deviation was computed to determine the 
closeness of fit to the measured points of the various regions. 
The closeness of fit is an indication of the error involved. The 
least standard deviation was found to be when the data was 
utilized along as much of the forward curve as possible. A 
comparison using the standard deviation from data received in 
different regions along the curve reveals this fact. Using the 
points corresponding to the large values of the measured voltage, 
or far out on the curve, the standard deviation is three times 
greater than the standard deviation found using the complete 
curve. The deviation found from the data of small values of 
voltage, or closest to the origin, is also much larger than the 
least standard deviation. These calculations are included in the 
Appendix using the 1N540 diode. 

The reverse saturation current determined from the forward 
characteristics was) used to help find Rj from the reverse 


characteristics by setting I, a constant an the derivati@n .«§ This 
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method reveals that this value of it has very little effect on 
the value of R, obtained. A simplified equation for determining 
Rj can be made, 
i a V5/ (14-5) (4.24) 

where V5 and I; are measurements from the reverse characteristics. 

The parameters for two diodes, 1N540 and 1N277 are calculated 
in the Appendix. The standard deviation for each diode using the 
parameters determined is also given. The standard deviation was 
found using the measured voltage V5 and the calculated voltage VC; 
from the measured current I,. Equation 4.9 is solved placing the 
estimations of the parameters in the equation for the constants 
and using the various values of measured current. The value 


obtained is called VC; - The standard deviation is found by the 


following relationship 


standard deviation (ero) 





where N is the total number of data points used. This value 

of the standard deviation is used as an error indicator. The 
method with the largest standard deviation has the largest error 
using the same data points. For the two diodes used, it was 
found that the least-square-error method had a smaller standard 
deviation than the normalized least-square-error method, for 
both diodes. The difference is quite small and is not a good 
indicator of which method is best. From using these programs, 


it appears that the normalized least-square-error method has an 
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advantage in calculations when the measurements are at relatively 
great distances from the actual desired curve. Either method 
estimates parameters that give a curve that has less deviation 
than that expected due to variations from diode to diode due to 


manufacturing procedures and techniques. 


B. DETERMINATION OF PARAMETERS FOR JUNCTION TRANSISTORS 

The model for the junction transistor was presented in 
Figure 2. The parameters to be determined are Ira, Ica, Me, Me, 
Rpp> Rcc> Rep Rp and Rc. In this analysis, the emitter-base and 
the collector-base junctions are considered to be two separate 
diodes for these calculations. For the emitter-base p-n junction 
parameters, measured data are made with the collector terminal 
open circuited and for the collector-base p-n junction parameters, 
measured data are made with the emitter terminal open circuited. 

Considering the emitter-base junction as a diode, calculations 
are made on the forward measured data to determine the quantities 
Tes, Me, and Rsf, the same as Is, M and Rs were determined for the 
discrete diode above. The quantity Rop is the combined emitter- 
base body, lead and contact resistance. The emitter reverse 
leakage resistance (Rp) is determined using the reverse character- 
istics, the same as for the diode. 

For the collector-base junction considered as a diode, 
calculations are made on the forward characteristics to determine 
the quantities Ic¢g, Mc, and Roc, the same as was computed for the 


diode. The quantity Roc is the combined collector-base body, 
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lead and contact resistance. The collector reverse leakage 
resistance (R,) ‘is determined using the reverse characteristics, 
the same as for the diode. 

To complete the model parameters for the junction transistor 
an assumption is made that the emitter bulk resistance is zero. 
The quantity Rep is the emitter bulk resistance. Using this 


assumption, Rpp is calculated from 


Rpg = (pp + DRop (4.26) 


The collector bulk resistance is determined to be 


the quantities HEE and HEE; are the common-emitter DC normal 
and inverted current gains respectively. 

The shift of ope can be evaluated from the specification 
sheets available on the individual transistors. Specification 
Sheets normally give curves ot Hee versus Ip- An average value 
is obtained by taking n values of HEE from the curve at equally 
spaced values of collector current over the range and then 
dividing by n. This average value of HFE, can then be used in 
further calculations of parameters. The parameter Hpp; can be 
measured directly or assumed to be values as shown in Table I. 

The common-base &'s can be calculated using equations 2.10 
anid e2e 11. 

This procedure of determining the model parameters for the 
Junction transistor are currently used by North American 


Rockwell Corp., among others. 
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TABLE I 


Typical Values of Hees 
(Reprinted from the Motorola Switching Transistor Handbook.) 


Transistor Type Typical HrEe; 
Planar, silicon epitaxial WAZ 
Alloy, uniform base u 
Alloy, diffused base “a : 
Mesa, diffused base 0.4 
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C. TYPICAL VALUES FOR MODEL PARAMETERS 

For the p-n junction diode the values of M, R. and ee Ge 
usually determined to be in the following ranges. The value of 
M is between 1.0 amd 3.0. the value of Ry is between 107? and 
10-1 kilohms, the value of I. silicon is from 10712 to 107°mA and 
for germanium is from 1078 to 10-*mA. The quantity Ry has a value 
between one and ten megohms. 

The junction transistor parameters are usually contained in 
the following ranges: The value of Mp and Mc is between 1.0 


= and 


and 3.0, the value of Rpp and Rec is between 10° 
107° kilohms, and Rpp is usually between iio and 1071 kilohms. 
Rpp has been assumed to be zero in this analysis. The usual 
range of Ipc is 10-14 +o 10-®mA for silicon and 1078 to 107°maA 
for germanium. The value of I¢c is usually an order of magnitude 
larger than Ipc . The quantities and Rpare usually larger 
than one megohm. 

These are only typical values, but if a parameter is 
substantially outside the usual range, the input data and 
calculations should be checked for error. These values were 
presented by Design Automation, Inc. Eas 

The parameter values for the two diodes calculated in the 
Appendix conform to the typical values for the model as is shown 
me caole 1. 

In Table III a comparison of the values calculated using the 
least-square-error methods to the values calculated following 
the procedures given by Design Automation, Inc. [i]. There is 
very little difference in the parameter values derived by either 


method. The standard deviation for the least-square-error methods 
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TABLE II 


INS4HO 
Usual Normalized 
deg fed 1.76 


1.65x10729 = = 1.86x19-10 
Te 0.126 


Q.006 0.0062 
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Calculated Parameters and Standard Deviations 


1N277 
Usual Normalized 
al. G7 1.10 
2.65x10719 = 3..72x19710 
82.83 67.9 
0.0025 0.0027 


TABLE III 


A Comparison of Parameters Derived Using Two Methods 


Design Automation 
AH di 


1.5x10719 
ieee 


0.0076 


Least Square 
ee 


1.65x10719 
0.131 


0.006 
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Normalized Least Square 
LPs) 


1.86x10-10 
Tews 


0.0062 


is a little less than that for the method of Design Automation, 
Inc. This indicates that the least-square-error methods are as 


well suited to these models as other procedures. 


D. COLLECTION OF DATA POINTS 

The measured values of current and voltage used in the 
A matrix and the b vector for the least-square-error methods need 
to be evenly distributed from the minimum value of voltage (at 
least 0.1 volt) to a value of voltage that corresponds to maximum 
power. If these values are too close together or on only one 
portion of the curve, the derived parameters turn out to be 
negative values or values that do not correspond to the usual 
range. Most of this type of error is due to the inaccuracy of 
the measurements of the voltage and current. A typical derivation 
is presented in the Appendix. After the parameters have been 
initially estimated, the additional points are utilized in the 
PECUNG Ive. Pellet MaMeKE Tie. The additional points may be taken 
anywhere between the minimum and the maximum values of voltage. 
The least standard deviation is found utilizing as much of the 
forward characteristics as possible. 

For the reverse characteristics, fewer points are needed and 
they may be spread more than the forward characteristics. The 
reverse leakage resitance is a large value and if equation 4.24 
is used only one set of points is needed in the calculation. 

This calculation need not be too accurate since the reverse 
leakage resitance is of little importance. Values of current 


and voltage between -0.2 volts and the breakdown voltage are 


acceptable. 


U2 


V. CONCLUSION 


"A good model has been chosen since it has been shown that 
these models accurately represent electrical performance of the 


device over a very wide range and the models are computationally 


efficient in regards to the computer circuit-analysis programs 


' 
§ 


presently in use. It is found that the error by ce ictal is 
considerably less than the variation in the characteristics due 
to manufacturing techniques. The region ieee to measure the 
data for the parameters for the forward eee rer ig | 
represented as the points lying between 0.1 volts and the maximum 
power rating of the device. The measured data best represents 
the device when the points cover as much of the curve between the 
above mentioned points as possible. Considering the two methods 
of parameter approximation, both methods give results that are 
well within the accuracy expected. The normalized least-square- 
error solution does more to eliminate the influence of points 
located farthest from the desired curve. 

The standard deviation, which is a measure of the error in 
parameter estimation, reflects that these models give an accurate 
representation of the actual performance of the semiconductor 
devices. The standard deviation from the least-square-error 
methods is a few millivolts and most measurements could be this 
inaccurate. As seen on the forward characteristics presented in 
Figures 6, 7 and 8, the difference between the actual measurements 


and computed values is almost indistinguishable. 
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Parameters were not estimated for a transistor, but the 
method of determining these parameters is straight-forward and 
should not present any serious difficulty. Semiconductor devices 
with more p-n junctions than the transistor can be modelee 
similarly by dealing with each junction separately as a diode. 

A comparison of values of parameters from the least-square 
estimates with other methods of parameter estimation reveals that 
the least-square-error methods produce parameter values essential- 
ly the same as other procedures. The error for the least-square- 


error methods is slightly less than other procedures. 
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‘ APPENDIX 


Examples of Parameter Estimation 


This section is used to give two examples of the use of the 
procedures presented in this thesis. Two diodes, one the LN540 
and the eee the 1N277, were used fn ths analysis as _ the two 
examples. Measured data of current and Heueer were taken in the 
forward and reverse Beans of ane two devices to — eee to 
aren ne Reeu device parameters. 

The following is the forward data collected for the two 


Peyiiee s with the voltage given in volts and the current in 


. “4 ’ 


Sees iekaperes 


IN540 LN27/7 
Voltage Current Voltage Current 
0.4 0.002 O225 0.002 
0.48 0.006 O..27 0.004 
0.5 0.0092 0.28 0.006 
0.52 0.014 O.2285 0.008 
0.54 Oae2 1 0.29 0.0L 
0.56 0.033 0.3 0.014 
Oegod § 40555 0.31 0.02 
0.6 0.086 ®.32 0.025 
0.64 e215 0.34 0.049 
0.68 0.483 0.36 0.085 
O.%2 ees 0.38 O.1L4 
0.76 3.18 0.4 0.222 
0.8 9.15 O.42 0.37 
0.84 oe) 
0.88 37.8 
0.92. 76 3D 
0.96 W550 
ibys ieee 270.0 


The reverse data for the devices are presented with both the 
current and the voltage as negative quantities and the current is 


in milliamperes and the voltage in volts. 


4S 


IN540 1N277 


Voltage Current Voltage Current 
9.0 0.01 2.65 0.004 
1 see O02 28 0.008 
Lee'0 0.014 10.0 0.0125 
5.5 0.016 ile 5 0.014 - 
mo. 5 0.20 20.0 0.024 
22 6 0.24 30.0 0.0355 
26.5 0.0265 36.0 0.0425 


Using seven selected points, an A matrix is assembled to 
start the recursive operation with the forward characteristics. 
These seven selected points are spaced over as much of the curve 
as possible so as to take into consideration the influence of 
Aire nente parameters on various regions of the curve. Errors are 
sometimes encountered when two points are used that lie too close 
to each other. 


The seven points chosen for the A matrixsand Devector sane 


IN540 1N277 
Voltage Current Voltage Current 
0”. 58 0.0555 0.25 0.002 
0.6 0.086 0.27 0.004 
0.68 0.483 0.285 0.008 
0.72 1 25 0.3 0.014 
On 86 37 8 0.34 0.049 
0.96 T5520 0.38 0.14 
1.0 270.0 0.42 0.37 


where the voltage is in volts and the current is in milliamperes. 
The natural logarithms of the currents were determined and the 
information was then run on the computer programs to start the 
recursive procedures. The solution of the parameters using only 
the first seven points for the 1N540 are as follows for the 


least-square-error methods described above. 


Yo 


Least Square Normalized Least Square 


A3 1.0249 1.0233 

Ap 0.0453 0.0451 

Rs 0.1268 0.1324 

Ay 6 .656x107 7. 1373x000" 
ie 1.502x10-10 1.4019x10-10 
M 1.743 1.734 


The solutions for the points for 1N277 are 


Least Square Normalized Least Square 
A3 0.6097 0.6162 
Ao 0.0275 peoZsn 
Rs 78.2493 71.3765 
Al 4.1626x109 3.3997x109 
he 2 tig2socheeee 2.941x10-10 
M 1.0587 1.0798 


The rest of the measured points are used in the recursive 


equations and result in the final forward parameters as 





1NS40 
Least Square Normalized Least Square 
A3 1.0246 1.0261 
Ao 0.0455 0.0458 
ne 0.1308 0.1257 
Al 6 .0536x102 5.3908x109 
ie 1.6519x107-10 1.855x10-10 
M 1 7s 17 ol 
1N277 
Least Square Normalized Least Square 
A3 0.6116 0.6236 
Ao 0.0277 0.0287 
Re 82.8309 67.8968 
Al 3.7768x102 2.6885x102 
ite A ease Lb 3.7195x10~-10 
M 1.0666 ote 


The seven points from the reverse characteristics were only 
evaluated using the normalized least-square-error method since 


the only parameter attained from the reverse characteristics is 
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the reverse leakage resistance R,;. The value of R, is of little 
importance and any value in the megohm range will be adequate,. 
The calculated values of R, are 

1N540 N27 


R 9.65x10”  8.69x107 


ut 
Values of these parameters were calculated with the procedures 
of Design Automation, Inc. [iJusing the data as listed above. The 


following-values were obtained. 


-LN540 
M ims 
Roose 


Fs 1.5x10719 


It is evident that these values are quite similar to those of the 
least-square-error methods. 
The standard deviation using equation 4.25 was calculated for 


the three sets of parameters calculated by the above methods.’ 


es 


Each standard deviation calculated is listed below. 
len eee 
Standard Deviations 


1Nn540 
Design Automation Least Square Normalized Least Square 
je 2? S07 0.006 0.0062 
1N277 
0.0025 0 G87 


In checking the various regions of the curve to be used the 
standard deviation was calculated for the parameters derived 
from data far out on the curve. The parameters calculated were 


M equals 1.872, Is equals 5.6@0xl0-19 and Rq equals 0.101. The 
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standard deviation for these parameters is 0.0184. This 
indicates that the best parameters are derived from data covering 
as much of the curve as possible. 

The standard deviation is a measure of the error between the 
calculated and the actual measurements. It is seen that the 
deviation is usually only a few milliamperes. 

Figures 6, 7, 8, 9, 10 and ll are plots of the forward 
characteristics from the measured data and also of the calculated 
parameters for diodes. From these curves it is clear that the 
solution obtained by the least-square-error methods is very 


accurate. 
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Plot of Measured Points for the 1N540 Diode 
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Plot of Calculated Points Using the Least-Square 
Parameters for the 1NSHO Diode 
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Plot of Calculated Points Using the Normalized Least-Square 
Parameters for the 1NS540 Diode 
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FIGURE 9 


Plot of the Measured Points for the 1N277 Diode 
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FIGURE 10 


Plot of Calculated Points Using the Least-Square 
Parameters for the 1N277 Diode 
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FIGURE 11 


Plot of the Calculated Points Using the Normalized 
Least-Square Parameters for the 1N277 Diode 
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CGC GleG) CO) 


COMPUTER PROGRAMS 


THIS IS THE MAIN PROGRAM FOR THE NORMALIZED LEAST a 
ERROR SOLUTION. 

A IS THE M BY N INPUT MATRIX OF COEFFICIENTS. 

B IS THE M BY 1 INPUT VECTOR OF CONSTANTS. 


RS, 


SI, AND ETA ARE THE OUTPUT PARAMETERS WHERE RS 


IS BODY RESISTANCE, SI IS THE REVERSE SATURATION CURRENT 


AND 


Ww 


118 


ee 


114 


25 


30 
32 


ETA IS THE VALUE FOR THE EMISSION CONSTANT M.. 
IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION A(10,10), B(10, 10), C(10,10), D(10,10) 
DIMENSION W(10,10), T(10), A2(10,10), E(10,10) 
DIMENSION A3(10,10), F(10,10), G(10,10), H(10) 
READ(5,4) M,N 

READ (5,5) ((A(I, od) , 11M) .d=1 No 
READ(5,5) (B(I,1) ,I=1 M) 

FORMAT (2110) 

FORMAT (7F10.0) 

DATA T/10*0.0/ 

DO 118 I=1,M 

DO 118 J=1,N 

T (1) =T (1) +A (1, J) **2 

CONTINUE 

DO 122 I=1,M 

W(I,1I) =1.0/DSQRT (T (I) ) 

CONTINUE 

DO 114 I=1,M 

B(I,1)=B(I,1) *W(1I,1) 

DO 114 J=1,N 

A(I,J) =A(1I,J) *W(1,1) 

CONTINUE 

DO 23 I=1,M 

nO. 23 J=1L.N 

A2(J,1)=A(1I.,J) 

CONTINUE 

CALL MATMUL(A2,N,M,A,M,N,C,M4,N4,L) 
CALL GAUSS3(N,EPS,C,A3,KER, 10) 

TP (KER. EO?) Gi. Tom 

CALL MATMUL (A3 ,M4,N4,A2,N,M,D,M5,N5,L2) 
CALL MATMUL(D,M5,N5,B,M,1,G,M8,N8,L3) 
PRINT 8, (G(I,1) ,I=1,N) 
FORMAT (10X, 3D14. 6) 

DO 303 J=1,10 

READ (33 0) o0( HiT) El) 

READ (5, 32) 

FORMAT (3F 10.0) 

FORMAT (1F 10.0) 

CALL WRLSF (A3,H,G,Z,N) 

PRINT 8, (G(I,1) ,I=1,N) 

A4=G (1,1) 
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ft) GG Gao 


33 


305 


10 


300 


AO=G (2,1) 

RS=G (3,1) 

A1=DEXP (A4/A0) 

PRINT 33, A4,A0,RS,A1 
FORMAT (LOX, 4D14. 6) 
SI=1.07AL 
ETA=A0*38 .46 

PRINT 8, SI,ETA 
CONTINUE 

GO TO 300 

WRITE (6,12) 

FORMAT (LOX, SHSINGULAR) 
STOP 

END 


THIS IS THE MAIN PROGRAM FOR THE LEAST SQUARE ERROR 
SOLUTION. 

A IS THE M BY N INPUT MATRIX OF COEFFICIENTS. 

B IS THE M BY 1 INPUT VECTOR OF CONSTANTS. 


RS, 


SI, AND ETA ARE THE OUTPUT PARAMETERS WHERE RS 


IS BODY RESISTANCE, SI IS THE REVERSE SATURATION CURRENT 


AND 


Zo 


Si) 
a2 


ETA IS THE VALUE FOR THE EMISSION CONSTANT M. 
IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION A(10,10), B(10,10), C(10,10), D(10,10) 
DIMENSION A2(10,10), A3(10,10), G(10,10), H(10) 
READ(5,4) M,N 

Pee, eeucl ,d) Sl Me d-1 SN) 
READ(5,5) (B(1,1) ,I=1,™) 

FORMAT (2110) 

FORMAT (7F 10.0) 

DO 23 I=1,M 

DO 23 J=1,N 

Alhaji 00k 1) 

CONTINUE 

CALL MATMUL(A2,N,M,A,M,N,C,M4,N4,L) 
CALL GAUSS3 (N,EPS,C,A3,KER, 10) 

IF (KER.EQ.2) GO TO 10 

CALL MATMUL(A3,M4,N4,A2,N,M,D,M5,N5,12) 
CALL MATMUL(D,M5,N5,B,M,1,G,M8,N8, L3) 
PRINT 8, (G(I,1) ,I=1,N) 
FORMAT (10X, 3D14.6) 

DO 303 J=1,10 

READ(5,30) (H(I) ,I=1,N) 

READ(5,32) Z 

FORMAT (3F 10.0) 

FORMAT (1F 10.0) 

Chl, TRESE(ASH,G,Z,N) 

BEGNI a. (G(1.1),1=1,N} 

A4=G(1,1) 

AO=G (2,1) 

RS=G (3,1) 

A1=DEXP (A4/A0) 


a7. 


33 


303 


10 


300 


Jie 


52 


124 


26 


134 


128 


130 


PRINT 33, A4¥,A0,RS,AL 
FORMAT (10X, 4D14.. 6) 
SI=L.O/AL 
ETA=A0¥*38 .46 

PRINT 8, SI,ETA 
CONTINUE 

GO T@ 300 

WRITE (6,12) 
FORMAT (10X, 8HSINGULAR) 
STOP 

END 


SUBROUTINE WRLSF (A3,H,G,Z,M) 
IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION A3(10,1) ,H(1), Q(10) ,v2(10), G(10,1) 
DIMENSION Q2(10,10) ,03(10,10), G1(10) 
CALL MATVEC (A3,M,M,H,Q) 

CALL VECMUL(H,H,M,S) 

CALL VECMUL(H,0,M,V) 

V1=S+V 

DO 122 I=1,M 

V2 (L) =O'(G)4V1: 

DO 132 I-1,M 

Glia 1D) 

CONTINUE 

CALL VECMUL(H,G1,M,Z2) 

7,3=Z-Z2 

Do! 12e1=1,M 

V2 (I) =V2 (I) *23 

DO 126 I=1,M 

G1(I) =G1 (1) +V2 (1) 

DO 134 I=1,M 

a1.) Seag' 

CONTINUE 

CALL MULVEC(M,0,H,Q2) 

CALL MATMUL(Q2,M,M,A3,M,M,03,M,M,L9) 
Q4=1.0D0+V 

DO 128 I=1,M 

DO 128 J=1,M 

Q2(I,J) = Q3(I,J) “ou 

DO 130 I=1,M 

DO 130 J=1,M 

Ade) =A) Oe (1a) 

CONTINUE 

RETURN 

END 
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ee 


132 


124 


126 


134 


26 


lO 


SUBROUTINE TRLSF(A3,H,G,Z,™) 
IMPLICIT REAL*8 (A-H,0-Z) 

DIMENSION A3(10,1) ,H(1), Q(10), V2(10), G(10,1) 
DIMENSION Q2(10,10) , Q3(10,10), G1(10) 
CALL MATVEC (A3,M,M,H,0) 

CALL VECMUL(H,Q,M,V) 

V1=1.0+V 

DO 122 I=1,M 

V2 (I) =O (I) W1 

DO 132 I=1,M 

Gilly Ge) 

CONTINUE 

CALL VECMUL (H,G1,M,2Z2) 

Z3=9PZ2 

DO 124 I=1,M 

V2 (I) =V2 (1) #23 

DO 126 I=1,M 

G1(I) =G1 (1) +V2 (I) 

DO 134 I=1,M 

ele) Cts) 

CONTINUE 

CALL MULVEC(M,0,H,Q2) 

CALL MATMUL(Q2,M,MA3,M,M,03,M,M,L9) 
DO 128 I=1,M 

DO 128 J=1,M 

24a) Oona) 7 

DO 130 I=1,M 

DO 130 J-1,M 

A3 (I,J) =A3(I,J) -02(1,0) 

CONTINUE 

RETURN 

END 


SUBROUTINE VECMUL(X,Y,M, SUM) 
IMPLICIT REAL*8 (A-H,0-Z) 
DIMENSION X(1) ,Y(1) 

SUM = SUM-SUM 

DO 5S I=1,M 

SUM=SUM+X (I) *Y (I) 

CONTINUE 

RETURN 

END 
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SUBROUTINE MULVEC (M,X,Y, XY) 
IMPLICIT REAL*8 (A-H,0-Z) 
DIMENSION X(1) ,Y(1) ,XY(10,1) 
DO 6 I=1,M 
DO 5 J=1,M 

5 XY(I,J) =X(I) *Y(J) 

6 CONTINUE 
RETURN 
END 


SUBROUTINE MATVEC (A,M,N,V,P) 
IMPLICIT REAL*8 (A-H,0-Z) 
DIMENSION A(10,1) ,V(1) ,P(1) 
DO 6 I=1,M 
Pel y=Pa) -P()) 
DO 6 J=1,N 
P(I) =P(I) +A(I,J) *V(J) 

6 CONTINUE 
RETURN 
END 


SUBROUTINE MATMUL(A,MA,NA,B,MB,NB,C,MC,NC,L) 
IMPLICIT REAL8* (A-H,0-2Z) 
DIMENSION A (10,1) ,B(10,1) ,C(10,1) 
L=0 
IF (NA.NE.MB) RETURN 
L=1 
MC=MA 
NC=NB 
DO 12 J=1,NC 
DO 12 I=1,MC 
e080 
DO 12 K=1,NA 
12 C(I,J)=C(1,J) +A(1I,K) *B(K,J) 
RETURN 
END 


The other subroutine GAUSS3 is a subroutine in the Naval 
Postgraduate School Computer Facility subroutine library for 
calculating the inverse of a matrix if such an inverse exists. 

These programs are not the most computationally efficient but 
are adequate for this calculation. There are several improvements 


that can be made to these programs expecially in regard to 


storage required. 
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